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ABSTRACT 

■ Radial trends of stellar populations in galaxies provide a valuable tool to understand 

the mechanisms of galaxy growth. In this paper, we present the first comprehensive 
analysis of optical-optical and optical-NIR colours, as a function of galaxy mass, out 
to the halo region (8R e ) of early- type galaxies (ETGs). We select a sample of 674 
massive ETGs (A/* > 3 x 10 10 Af Q ) from the SDSS-based SPIDER survey. By com- 
paring with a large range of population synthesis models, we derive robust constraints 
on the radial trends in age and metallicity. Metallicity is unambiguously found to 
decrease outwards, with a measurable steepening of the slope in the outer regions 
(Re < R < 8R ). The gradients in stellar age are found to be more sensitive to 
the models used, but in general, the outer regions of ETGs feature older populations 
compared to the cores. This trend is strongest for the most massive galaxies in our 
sample (M* ^ IO^Mq). Furthermore, when segregating with respect to large scale 
environment, the age gradient is more significant in ETGs residing in higher density 
regions. These results shed light on the processes leading from the formation of the 
central core to the growth of the stellar envelope of massive galaxies. The fact that 
the populations in the outer regions are older and more metal-poor than in the core 
suggests a process whereby the envelope of massive galaxies is made up of accreted 
C3 ' small satellites (i.e. minor mergers) whose stars were born during the first stages of 

galaxy formation. 

Key words: galaxies: fundamental parameters - galaxies: elliptical and lenticular, 
cD - galaxies: evolution - galaxies: formation - galaxies: stellar content - galaxies: 
groups: general 
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1 INTRODUCTION 

The link between chemical and dynamical evolution in 
early-type galaxies (hereafter ETGs) has been studied for 
quite a long time through systematic measurements of 
colour gradients, which may lead us to understanding how 
the variations in the age and metallicity of the underly- 
ing stellar populations evolve as st ar formation proceeds 



through the history of the system dd e Vaucouleurs 1961; 
Davies et al.ll 19871; [Peletier et al.l | l99d ; lFerreras et al.ll2u05a ; 
La Barbera fc de Carvalhq 20091 ') . The theoretical models 



suggest that metallicity and age gradients arise naturally 
from the processes leading to galaxy formation, li ke in the 
"in-s itu" collapse of a proto-galactic gas cloud (|Carlberd 
I1984T ). The first rendition of th is formation scena r io, also 
known as monolithic collapse jEggen et all Il962l ; ILarsonl 
1 19741 . Il975l ). predicted steep metallicity gradients, while, a 
revised monolithic model for the form ation of ETGs led 
to sh allower metallicity gra dients (e.g. Pipino et al.l 1 20081 . 
l20ld V Simulations (see e.g. lKobavashill2004 l indicate that 
galaxy mergers will weaken the metallicity gradient, and 
lead to significant levels of star formation in the galaxy cen- 
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tre. More recently, cosmological simulations have brought to 
our attention the importance of cold accretion, as a viable 
way to form galaxies " in-situ" , resemblin g the monolithic 
formation scenario (e.g. iKeres et al.ll2005l ). It is still a mat- 
ter of debate whether the metallicity gradient we measure 
today is a product of an old stellar population - reflecting 
mainly the initial conditions when the system collapsed - or 
results from more recent star formation episodes. Therefore, 
systematic measurements of colour gradients are of crucial 
importance to distinguish among different formation models 
of early-type systems. A proper constraint on the gradients 
of the properties of the underlying stellar populations will 
help understand the mass assembly process in ETGs. 

In the past, spectroscopic indices were used to mea- 
sure metallicity in ETGs, revealing the existence of ra- 
dial gradients rang i ng from —0-1 to —0.3 dex per decad e 
jCarollo et al.ll 19931 ; iDavies et alll993l ; iMehlert et ai]|2003h . 
The number of ETGs for which the radial dependence of 
age and metallicity extends beyond one effective radius is 
small for today standards, albeit growing. This is due to the 
difficulty in obtaining spectroscopic measurements at sur- 
face brightness below the background sky level. However, 
the availability of multi-waveband samples for a vast num- 
ber of sources remedy this situation, allowing us to esti- 
mate colour gradients up to several e f fective radii at low 
to m oderate redshifts l|Wu et al.1 120051 ; iTal fc van Dokkuml 
l201ll ), and op en the possibil i ty of studying radial gradi- 
ents at z>l l|Gargiulo et al.1 l201ll ; iGuo et al.l l201ll ). An 
important caveat is that colours have to be modelled by 
a stellar population synthesis code. This approach has 
been used in the past for small sam ples of ETGs (e.g. 
iPeletier et "all Il990l : ISaglia et al.ll2000l ). while in the recent 
years, the study of large gala xy samples has mostly focused 



on the optical reg i me (e.g. iRoche. Bernardi Hvdej l20ld : 



iTortora et al]|2010l : iGonzalez- Perez 1 et al.ll201lir where one 



has t o contend with t he age-metallici t y-ext inction degener- 
ISilva fc Elstonl (| 19941) . and more re- 



acy JWorthevlll994l) 

;lv iRoediger et al.l |201ll ). 



cently Roedig er et al.l (|2011f ). have found evidence for sig- 
nific ant age gradients in ellipticals, while ind ependent works 
fe.g. lHinklev fc Imll200ll ; lMehlert et ai]|2003l ), found negligi- 
ble age gradients, showing that the interpretation of colour 
gradients is a difficult task. 

Another aspect of studying colour gradients in ETGs 
is to understand how this quantity relates to other observ- 
ables and relationships so that we can build up a consistent 
galaxy formation and evolution framework. For instance, al- 
though ETGs obey specific scaling relations (e.g. Funda- 
mental Plane), they may well represent different families ac- 
cording to the way mass was assembled during their history 
and the characterization of how mergers may have estab- 
lished the mass as sembly is of unquestionable interest (e.g. 
I Greene et al.|[2012T ). As a first approximation, we would ex- 
pect mergers to significantly affect the well-known scaling re- 
lations involving central velocity dispersion and stellar pop- 
ulation properties, unless mass is accreted at very large radii. 
Theoretical and observational work has shown that the pres- 
ence of tidal debris in the outskirts of ETGs may provide 
an essential piece of evi dence to distinguish a mong different 
assembly scenarios (e.g. iDuc fc Renaudll201ll ). So far, most 
of the stellar population work was carried out within the 
central regions and mainly restricted to the optical regime. 
Extending the analysis to the outer regions with a combi- 



nation of optical and near-infrared data is a key factor to 
better comprehend how galaxies form and assemble their 
baryons. 

In this paper we combine optical (Sloan 
Digita l Sky Survey-Data Rel ease 6, SDSS- 
DR6; lAdelman-McCarthv et all [2008) and near-infrared 
(UKIRT Infra red Deep Sky Sur vey-Data Release 4, 
UKIDSS-DR4; lLawrence et al.1120071 ) data to study colour 
gradients in a sample of 674 ETGs located in different envi- 
ronments. The data were analysed with a dedicated pipelin e 
that uses the package 2DPHOT (|La Barbera et al.1 [2008). 
The high quality image of both optical and NIR surveys 
enabled us to determine reliable colour gradients out to 
eight effective radii, as shown below. All systematic effects 
were carefully taken into account, including background 
subtraction; stacking of the colour profiles and the effect of 
the tails of the PSFs in the redder bands, all of which are 
crucial for an accurate measurement of the colours in the 
outer regions. The results presented here confirm previous 
findings about the metallicity gradients, and show that age 
gradients may not be as negligible as some papers have 
reported in the past. 

The layout of the paper is as follows. In Sec. [2] we present 
our sample of ETGs. Sec.[3]deals with the methodology fol- 
lowed to determine the colour profiles, and discusses possi- 
ble sources of systematics. In Sec. [4] we present the median- 
stacked colour profiles of ETGs, from g — r through g — K, 
out to a maximum galactocentric distance of 8R e , for dif- 
ferent galaxy mass bins. Sec. [5] describes the fitting of the 
observed profiles with theoretical models based on different, 
state-of-the-art, stellar population synthesis codes. In Sec.[S] 
we show the results of the stellar population fits, i.e. how 
age and metallicity are found to vary from the centre to 
the external regions of ETGs. We also discuss the contribu- 
tion of internal reddening. Sec.Qshows how the inferred age 
and metallicity profiles are affected by galaxy environment, 
while in Sec. [S] we summarize the main findings of this paper. 
Throughout the paper, we adopt a standard ACDM cosmol- 
ogy with Ho^Skms -1 Mpc~\ f2 m = 0.3, and Qa = 0.7. 



2 THE SAMPLE 

The present study is base d on a set of 674 ETGs, extracted 
from the SPIDER survey (|La Barbera et al . 2010a, hereafter 
Paper I). This sample selects those galaxies with the best 
quality of the parameters that describe their surface bright- 
ness distributions. The SPIDER sample consists of 39, 993 
luminous ETGs, in the redshift range of 0.05 to 0.095, with 
M r < — 20, where M r is the k-corrected SDSS Petrosian mag- 
nitude in r-band. All galaxies have griz photometry and 
spectroscopy from SDSS-DR6, while 5, 080 ETGs also have 
photometry in the YJHK wavebands from the UKIDSS- 
Large Area Survey (see Paper I). In all wavebands, the 
structural parameters, i.e. the effective radius, R c , mean 
surface brightness within that radius, (/Lt) e , and Sersic in- 
dex, n, have bee n homogeneously measu red using the soft- 
ware 2DPHOT |La Barbera et al.1 120081 ). by fitting galaxy 
images with two-dimensional seeing-convolved Sersic mod- 
els. Total magnitudes are computed from R c and (fi) e in each 
filter. Stellar masses are derive d by fitting synthetic ste l- 
lar population (SP) models from lBruzual fc Charlotl (2003), 
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Figure 1. Comparison of the median surface brightness profiles of ETGs in different bins of stellar mass, M* (increasing from left to 
right), and axis ratio, b/a (decreasing from top to bottom). For each panel (i.e. a given bin of M* and b/a), black curves plot the r-band 
surface brightness profiles of all galaxies in that bin. The profiles are shown as a function of the normalized galactoccntric distance R/R c , 
where R is the r-band effective radius. For each panel, the solid red curve is obtained by median-stacking all the single profiles, with 
dotted curves marking the ±1<t scatter around the median profiles. For each bin of M*, the median profile that corresponds to highest 
b/a (top panel) is repeated in all panels (from top to bottom) as a coloured dashed curve (blue, magenta, and orange for left, middle, 
and right panels, respectively). 



with a variety of star formation histories and metallicities, 
to th e optical+NIR ph otometry, usi ng the softwar e LeP - 
hare (|llbert et alj |2006), assuming a ICardelli etlT] dl989h 
extinction law and Chabrier IMF (see ISwindle et all 20 111 
hereafter Paper V, for details on the estimate of M*). 

To perform an accurate stacking of the surface bright- 
ness profiles (seeing de-convolved) as a function of galacto- 
centric distance (see Sec. , we select ETGs according to 
the following criteria. 

- We select only galaxies whose light distribution is well 
fitted by a Sersic model, i.e. with x 2 < 2 in all wave- 
bands, where \ 2 ls the rms of residuals between the galaxy 
image in a given band and the corresponding best-fitting 
two-dimensional Sersic model. We also remove galaxies with 
large uncertainties (> 0.5 dex) in logR c , leading to a sam- 



ple of 4, 546 (out of 5,080) ETGs. Th is sample is the same 
one used in lLa Barbera et all l)2010d . hereafter Paper IV) , 
where the inner (< 1 R c ) colour gradients were analyzed, 
along with their dependence on mass and stellar population 
properties (i.e. age, metallicity, and [a/Fe]). 

- We select galaxies in a narrow redshift range, from 0.05 
to 0.07, which corresponds to those with the best S/N ra- 
tio, minimizing the effect of k- and evolutionary corrections 
varying within the sample and as a function of galactocentric 
distance. This results in 1,255 ETGs. 

- Galaxies with stellar mass less than M* = 3 x 10 10 M© 
are removed, leading to 1,043 ETGs. The reason for this 
selection is the following. At the upper redshift limit of 
z — 0.095, the magnitude limit of the SPIDER sample 
(M r ~ —20; see above), corresponds approximately to the 
magnitude limit below which SDSS spectroscopy becomes 
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incomplete (i.e. a Petrosian magnitude of m r ~ 17.8). This 
makes the SPIDER sample incomplete at a stellar mass of 
M* ~ 10 10 M Q (corresponding to M r ~ -20). The selec- 
tion limit of Mi, = 3 x 10 10 A/q ensures that the sample of 
ETGs used for the present analysis is 90% complete with 
respect to stellar mass, this completeness level being esti- 
mated from the M 4 versus r-band Petrosian magnitude plot 
(see, e.g., sec. 3.2 of Paper IV). 

- Finally we remove galaxies with an axis ratio, b/a ^ 0.65. 
The reason for this selection is illustrated in Fig[T] that com- 
pares the r-band median-stacked surface brightness profiles 
(solid red curves) of the 1,043 ETGs at 0.05 z 0.07 
(see above) for three different bins of b/a (from top to bot- 
tom) and (from left to right). For a given bin of M*, the 
median profile corresponding to the highest b/a is shown 
in all the panels as a dashed curve. At all masses, the pro- 
file does not change significantly (within a few percent) for 
b/a Si 0.65, justifying ou r cut based on b/a. As argued by 



iHvde fc Bernard! (|2009h . ETGs with b/a < 0.6 may indeed 
represent a different population of objects (e.g. rotationally 
supported systems). 

Our final sample of high-quality data comprises 674 
ETGs, and we split this sample in three bins of stellar 
mass, keeping the number of galaxies per bin fixed. Here- 
after, we refer to galaxies in the three bins as low- (3.00 < 
M*[10 10 A'/ Q ] < 5.29), intermediate- (5.29 < Af*[10 10 M@] < 
9.42), and high- (9.42 < Af*[10 10 Af Q ] < 70.71) mass, re- 
spectively. In each bin of galaxy mass, we also split galaxies 
according to the environment where they reside. To this ef- 
fect, we use the friends-of- f riends (FoF) group catalogue, as 
described in lBerlind et al. (2006). We use an updated cata- 
logue b ased on SDSS-DR7 - rather than DR3, used for the 
original iBerlind et al. I l|2006l ) dataset. The catalogue, that 
comp rises 10, 124 syst ems, is subject to a virial analysis 
fsee lLopes et al.ll2009al . and references therein), resulting in 
a compil ation of 8, 083 groups wi th well measured proper- 
ties (see lLa Barbera et aL I l2010bl . hereafter Paper III, for 
details). ETGs in the SPIDER sample are then classified 
as either (i) group galaxies, i.e. those having group mem- 
bership according to the virial analysis (~ 46%); (ii) field 
galaxies (~ 33%), i.e. objects located more than five virial 
radii away from any group initially detected by the FoF 
algorithm; and (iii) unclassified (~ 21%), i.e. objects in nei- 
ther of the other two classes. Out of the selected sample 
of 674 ETGs analyzed in the present work, 336 (205) ob- 
jects are group (field) galaxies. The group galaxies reside 
in systems with average velocity dispersion (virial mass) of 
- 240km s" 1 (~ 8 ■ 1O 13 M ). These subsamples are used to 
investigate the dependence of stellar population profiles in 
ETGs on environment (Sec. [7J| . 



3 DERIVATION OF COLOUR PROFILES 
3.1 Methodology 

For a given waveband X, with X — grizY J H K , the surface 
brightness profile of each ETG, fix(R), is modelled by the 
Sersic law: 



a x (R) = mo,x + 1.0857 • b„,x 



R 



Re 



1/»X 



where no.x is the central surface brightness of the Sersic 
model in the passband X; nx is the corresponding Sersic 
(shape) parameter; and fc„,jf ~ 2nx ~ 1/3 is a con- 
stant defined so that R e ,x is the half-light radius of the 
model {Caon~eFld][l993). The radius, R, is the circularized 
galactocentric distance . For each galaxy we construct the 
colour profile, in the form g — X (with X — rizYJHK), by 
computing the difference of the corresponding Sersic profiles: 



g—X (i?) = tto, 9 - t*o,x + 



+ 1.0857 • 



R 

Re.c 



l/r, 



bn,X 



R 



R 



e,X 



l/r, 



(2) 



All the parameters ^o,x , Re,x, and nx, are obtained by fit- 
ting the galaxy image in waveband X with a PSF-convolved 
Sersic model, as detailed in Paper I. For a given galaxy mass 
bin, a stacked colour profile is then computed. To this effect, 
for a given galaxy mass bin, the galactocentric radii of each 
individual colour profile are normalized by the correspond- 



ing R e 



and the median of all the so-normalized colour 



(1) 



profiles is computed. 



3.2 Systematics 

At the very faint surface brightness levels we want to explore 
(i.e. Si 25 mag arcsec -2 in the r band, at distances larger 
than a few effective radii, see Fig. [TJ, the colour profile of 
ETGs can be affected by several sources of systematics. We 
discuss here the role of such systematics, and how we take 
them into account in the analysis. 

Background subtraction. As described in Paper I, struc- 
tural parameters of ETGs are estimated by fitting galaxy 
images with a PSF-convolved Sersic model plus a con- 
stant value representing the background. An error in the 
background estimate can bias the R c ,x and nx parameters, 
changing the shape of the colour profiles (see Eq. [2}. To 
quantify the relevance of this effect, we have repeated the 
Sersic fitting in g and r bands, fixing the background to 
the value measured in the outermost regions of e ach postage 
stam p frame by applying the biweight statistics (Bee rs et al.l 
1990). We found that the stacked g — r colour profile does not 
change significantly when using a fixed-background fitting 
(with respect to the case where background is treated as a 
free fitting parameter), with a variation < 0.02 mag arcsec -2 
at the largest radii probed (R/R c = 8), i.e. smaller than the 
typical statistical error of the stacked profiles at these radii 
( Si 0.03 mag arcsec -2 ; see Tables[T] and[3J). In this work, 
we use structural parameters, and hence color profiles, de- 
rived by treating the background as a free fitting parameter 
(i.e. adopting the same procedure as in Paper I) . 

Parametric vs. non-parametric profiles. Our stacked 
colour profiles are obtained by a parametric approach, us- 
ing the best-fitting Sersic models in different wavebands to 
compute the colour profiles. This approach is well justi- 
fied by the fact that, on average, the light profiles of mas- 
sive ellipticals are well described b y a single Sersic mode l 
out to roughly eight effec tive radii l|Kormendv et al1l2009l : 
iTal fc van DokkurrJ l201ll ') . Moreover, the % 2 selection (see 
Sec.[5J allows us to exclude galaxies whose light profile is not 
reproduced by the Sersic model. These objects are mostly 
early-spiral contaminants (see Paper I for details). In App.lAl 
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we show the quality of our Sersic fits, by comparing the 
(stacked) best-fitting light profiles, derived on circular aper- 
tures, with those measured directly from the galaxy images, 
for all relevant wavebands and each mass bin. Parametric 
and non-parametric light profiles turn out to be fairly consis- 
tent, within uncertainties, supporting the robustness of our 
measurements out to large galacto-centric distances. A fur- 
ther test of our parametric approach is illustrated in Fig. [2l 
comparing, as example, the stacked g — r profiles obtained 
by the parametric and non-parametric approaches for high- 
mass ETGs. Our choice of g — r for the comparison, rather 
than other colour profiles, is justified by the fact that this 
colour is less sensitive to the effect from the wings of the PSF 
(see below). Also, because of the shorter wavelength baseline 
with respect to other colors, g — r provides a more critical 
test. For each galaxy, we measured the surface brightness 
profiles in the g and r bands directly from the images over 
the same concentric ellipses in both bands, all ellipses having 
the same axis ratio from the r-band Sersic fit and circular- 
ized radii equally spaced by 0.5pixel (i.e. ~ 0.03 R e for high- 
mass ETGs). For each galaxy, we then rescaled the radii 
of the ellipses by a factor 1/R e , where R e , r is the r-band 
effective radius of the given galaxy. The rescaled profiles 
were median-combined in both g and r bands, to generate 
the non-parametric g — r stacked color profile. The stacked 
non-parametric profile was then averaged over radial bins 
logarithmically spaced by 0.15 dex, the profile uncertainty 
in each bin being computed from the standard deviation of 
colour values in the bin. Fig. [2] shows that parametric and 
non-parametric profiles differ significantly only in the inner 
part, at R/R e < 0.2, as expected, since the non-parametric 
profile is not corrected for seeing effects, while they agree 
quite well out to a galactocentric distance of ~ 3 — 4R e . We 
decide to continue our study using the parametric approach 
(i.e. Sersic fitting) as the corresponding uncertainties in the 
surface brightness photometry are significantly smaller at 
large R/Re, with respect to the case where a non-parametric 
approach is adopted - a crucial aspect for the present study. 

PSF red halos. The shape of the PSF at large radii can 
alter the colour profile of a galaxy. As regards SDSS data, 
the PSF exhibits prominent wings at distances ^ 10", be- 
ing m ore pronounced in t he i band than in the g and r 
filters (|Bergval et alJl201Ch . If neglected, the wings in the 
redder bands of the PSF cause a spurious re d colour excess 
(in g — i and r — i) in th e galaxy outskirts |de Jonell2008l ; 
iTal fc van Dokkumll201lh . Our colour profiles are based on 
PSF-corrected structural parameters. The PSF is obtained 
by fitting images of stars with a superposition of 2D Moffat 
functions, truncated at a maximum radius, rpsF, where the 
PSF model reaches a mean flux level of 1% of its central 
value. In the previous papers of the SPIDER survey, rpsF 
is typically ~ 7 pixels (i.e. ~ 2.8") in all wavebands. This 
issue is not relevant for those papers, as only the properties 
within 1-2 effective radii were considered. However, when 
studying the profiles out to 8R e , this choice will underes- 
timate the effect of the wings of the PSF. To account for 
this, we have reprocessed the griz images of all the 674 se- 
lected ETGs using a large PSF model, with rpsF = 60 pixels 
(~ 24"). Hereafter, we refer to the original (rpsF = 7 pixels) 
and the updated PSF models as the small and large PSFs, 
respectively. Fig.[3]compares the median-stacked g — r, g — i, 
and g — z, colour profiles of high-mass ETGs with respect 




R/R e (r band) 

Figure 2. Median colour profiles of high-mass ETGs as ob- 
tained by the parametric (dashed curve) and non-parametric 
(solid curve) approaches (see text for details). Dotted grey (black) 
curves mark the ±1 a error on the parametric (non-parametric) 
profile. The bump in the inner (R/R c < 0.2) part of the non- 
parametric profile is due to the variation of the PSF between the 
g and r bands (the latter having a a smaller FWHM, produc- 
ing the red bump). The profiles are in good agreement within 
the scatter out to the largest radii probed in the present study 
(~ 8R C ). Notice how the scatter in the non-parametric profile 
increases significantly at large galactocentric distances. 

to the size of the PSF model. Using a small PSF makes the 
galaxy colour redd er in the e xtern al regi ons, beyond ~ !R e , 
in agr eement with Ide Jongl (|2008l ) andlTal fc van Dokkuml 
l|2011l ). As expected fro m the waveband dep endence of the 
PSF halo (see fig. 6 of iBergval et alj l2010h . the effect is 
stronger in g — i (and g — z) than g — r. In the inner regions, 
at R/R < 1, the profiles (and hence the inner colour gradi- 
ents) are not affected by the choice of PSF model. For the 
analysis presented in this paper, we have obtained all opti- 
cal colour profiles using the large PSF. Notice that UKIDSS 
PSFs are not affected from any extended wing, as oppose to 
the i- and z-band (SDSS) PSFs. Also, the large wavelength 
baseline probed by optical-NIR colors makes them insensi- 
tive to any PSF red halo effect. For this reason, we have 
decided not to reprocess the NIR wavebands, and use the 
Sersic profiles from our previous works. 



4 COLOUR PROFILES VS. GALAXY MASS 

Fig. [4] plots the median colour profiles, g — X (with X — 
rizYJHK), obtained by median-stacking the single profiles 
of all galaxies in each of the three stellar mass bins. For 
each panel, the black dashed line represents the slope of 
the colour profile in the radial range from 0.1 to lR e , i.e. 
the average colour gradient, V 9 _x, as estimated in Paper 
IV. Error bars are la uncertainties on median colours, for 
intermediate-mass ETGs, at different fiducial galactocentric 
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Figure 4. Median colour profiles, in the form of g — X (with X = rizYJHK), for the sample of 674 ETGs, as a function of the 
normalized galactocentric distance R/R c . For each panel, the blue, magenta, and orange curves are the median profiles for galaxies in the 
same stellar mass bins, as labelled (see bottom-right of the Figure). The black dashed lines have slopes equal to the peak values of colour 
gradient distributions of SPIDER ETGs, as reported in tab. 2 of Paper IV. The intercept of the dashed lines has been arbitrarily chosen 
to match the inner (R < R e ) colour profile of intermediate mass galaxies (magenta curve). The six vertical error bars mark the positions 
of 0.1, 0.5, 1, 2, 4, 8 effective radii that we select to analyze the variation of colour gradients in ETGs from their centre to the outskirts. 
The size of the error bars is given by the ±lcr error on the median colour profile of intermediate mass galaxies (magenta colour). Notice 
the remarkable upturn in the g — i and g — z profiles at R ^ R c for intermediate- and high-mass galaxies. 



distances, i.e. R/R c = {0.1, 0.5, 1, 2, 4, 8}. The values of 0.1 
and 1 R e define the range usually ado pted to estimate in - 
ternal colour gradients of ETGs (e.g. IPeletier et allll990l ). 
while the distances {2, 4, 8} R e are chosen to analyze the 
behaviour of colour gradients in the outermost regions. We 
adopt a maximum radius of 8R e , so that the photometric 
errors are kept small enough ( ^ 0.04 mag) to allow for a 
meaningful stellar population analysis (Sec. 15.2] ). Moreover, 
as noticed above, a single Sersic model might not be able 
to describe th e light content of ETGs a t larger galactocen- 
tric distances (|Tal fc van Dokkum|[2(5ll1 ). The median g — X 
colours (with X = rizYJHK), computed at the fiducial 
galactocentric distances, are reported, together with their 
uncertainties, in Tables [1] [2] and [3] for low-, intermediate- 
, and high-mass ETGs, respectively. Notice that both the 
median g — X colours and colour uncertainties are obtained 
by stacking the parametric color profiles (Eq. [2} . 

Both optical-optical and optical-NIR colour profiles be- 
come bluer towards the galaxy outskirts, up to ~ 8 R c . Over- 
all, with the exception of g — r and g — J, the profiles are 
reasonably linear over the entire log-radial range. In g — i, 
low- and intermediate-mass galaxies exhibit some deviation 
from the linear trend at R/R c ^ 3, perhaps reflecting some 
inaccuracy in the i-band PSF halo modeling at large dis- 
tances. However, the deviation is comparable to the size of 
the error bars, making the effect unimportant for the stellar 
population analysis (Sec. I5.2p . The curvature in the g — r 
profiles is small, considering error bars, while the g — J pro- 
file shows a strong curvature, not seen in other optical-NIR 
colours (g — Y, g — H, and g — K). This might be explained 



by systematics in the J-band structural parameters, due to 
the different quality and resolution of the J-band photom- 
etry (see Paper I for more details). For this reason, we do 
not use g — J to constrain radial trends of stellar properties 
(Sec. [52) . 



This paper constitutes the first comprehensive analy- 
sis of optical-optical and optical-NIR colours, as a func- 
tion of galaxy mass, out to the halo region of ETGs 
(8 R e ). On the other hand, for what concerns optical data 
al one, we can compa r e our results with those obtained 
bv lTal fc van Dokkuml (|201ll . hereafter TvDll), whose anal- 
ysis was based on stacking a large sample of 42, 000 lumi- 
nous red galaxies from SDSS, at redshift z~ 0.34. Fig. [5] 
compares our g — r colour profile, for high-mass ETGs, with 
their r — i stacked profile. Notice that, at z~ 0.34, the r — i 
colour roughly corresponds to g — r in the rest-frame, al- 
lowing for a direct comparison to our colour profiles. The 
r — i profile has been extracted from fig. 8 of TvDll (left 
panel) and normalized to match our average g — r within 
R e . Our profile and that of TvDll are in very good agree- 
ment within an effective radius, while at large distances, the 
TvDll stacked data appear slightly bluer (by ~ 0.03 mag), 
but fully consistent with our results, considering error bars 
and the different selection criteria (e.g.. redshift range) of 
our and TdVll samples. 
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Table 1. Median colours of low-mass ETGs at different galactocentric distances. 



R/R c 


g-r 


9-i 


9-z 


g-Y 


g-J 


g-H 


g-K 


0.1 


0.871 ±0.009 


1.264 ±0.011 


1.588 ± 0.011 


2.382 ±0.015 


2.882 ±0.015 


3.588 ± 0.016 


3.957 ±0.017 


0.5 


0.829 ± 0.005 


1.201 ±0.006 


1.479 ± 0.008 


2.245 ±0.010 


2.741 ±0.013 


3.389 ±0.012 


3.732 ±0.014 


1.0 


0.805 ± 0.006 


1.186 ±0.007 


1.450 ± 0.009 


2.178 ±0.012 


2.665 ±0.015 


3.311 ± 0.014 


3.649 ±0.016 


2.0 


0.782 ±0.008 


1.166 ±0.009 


1.431 ± 0.012 


2.130 ±0.016 


2.579 ±0.026 


3.242 ± 0.019 


3.576 ±0.026 


1.0 


0.740 ±0.015 


1.151 ±0.017 


1.431 ± 0.022 


2.067 ±0.029 


2.427 ± 0.043 


3.178 ± 0.031 


3.475 ±0.040 


8.0 


0.695 ± 0.028 


1.157 ±0.033 


1.412 ± 0.040 


1.964 ±0.048 


2.263 ±0.072 


3.091 ± 0.054 


3.372 ± 0.062 



Table 2. Median colours of intermediate-mass ETGs at different galactocentric distances. 



R/R c 


g-r 


g-i 


g- z 


g-Y 


g-J 


g-H 


g-K 


0.1 


0.873 ± 0.008 


1.284 ± 0.009 


1.588 ± 0.009 


2.413 ±0.011 


2.936 ±0.012 


3.615 ± 0.014 


3.984 ± 0.014 


0.5 


0.838 ± 0.005 


1.224 ±0.007 


1.500 ± 0.009 


2.251 ±0.011 


2.758 ±0.013 


3.401 ± 0.014 


3.750 ±0.016 


1.0 


0.819 ±0.006 


1.192 ±0.008 


1.465 ± 0.010 


2.174 ±0.014 


2.651 ±0.019 


3.305 ± 0.016 


3.659 ±0.018 


2.0 


0.792 ±0.009 


1.175 ±0.012 


1.437 ±0.014 


2.086 ±0.019 


2.539 ±0.032 


3.207 ±0.027 


3.546 ± 0.027 


4.0 


0.768 ±0.019 


1.173 ±0.022 


1.401 ± 0.024 


2.006 ± 0.032 


2.391 ±0.053 


3.099 ± 0.042 


3.426 ±0.043 


8.0 


0.732 ±0.034 


1.182 ±0.039 


1.383 ± 0.039 


1.899 ±0.050 


2.206 ± 0.089 


2.981 ± 0.063 


3.285 ±0.067 
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Figure 3. Comparison of g — r, g — i, and g — z, colour pro- 
files, when using a small (rpsF = 2.8", grey curves) and large 
(''PSF = 24", black curves) PSF model (see text for details). 
Solid curves mark median trends, while dashed curves are the 
±lcr uncertainties on the median. 
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Figure 5. Comparison of our g — r stacked colour profile with 
the r — i stacked profile obtained for ETGs at redshift z~ 0.34 by 
TvDll. Notice that at z~ 0.34 the r — i colour corresponds ap- 
proximately to g — r in the rest frame. The TvDll profile matches 
well our colour profile in the inner galaxy region (R < R e ), while 
at larger radii their profile appears bluer (by ~ 0.03 mag), al- 
though fully consistent within the error bars, with our results. 

5 CONSTRAINING THE PROPERTIES OF 
THE UNDERLYING STELLAR 
POPULATIONS 

5.1 Stellar population models 

We compare the observed colours of ETGs with the predic- 
tions of three different versions of a wel l established stellar 
popul ation (SP) synthesis code, (a) the iBruzual fc Charlotl 
l|2003l. hereafter BC03) models , (b) the minor revision 
of these models introduced by ICharlot fc Bruzuall (2007, 
henceforth CB07), and (c) the major revision of this code 
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Table 3. Median colours of high-mass ETGs at different galactocentric distances. 



R/R c 


g-r 




g-z 


g-Y 


g-J 


g-H 


g-K 


0.1 


0.879 ± 0.005 


1.286 ±0.006 


1.610 ±0.007 


2.434 ±0.010 


2.965 ±0.012 


3.639 ± 0.013 


4.027 ±0.014 


0.5 


0.844 ± 0.005 


1.245 ±0.006 


1.529 ± 0.008 


2.275 ±0.011 


2.792 ±0.014 


3.435 ± 0.013 


3.789 ±0.015 


1.0 


0.822 ± 0.006 


1.225 ±0.007 


1.500 ±0.009 


2.200 ±0.014 


2.650 ±0.020 


3.355 ± 0.017 


3.699 ± 0.020 


2.0 


0.801 ±0.010 


1.205 ±0.011 


1.473 ± 0.015 


2.110 ±0.023 


2.492 ± 0.034 


3.267 ±0.027 


3.604 ±0.029 


1.0 


0.759 ±0.020 


1.193 ±0.022 


1.416 ± 0.026 


2.009 ± 0.041 


2.336 ±0.063 


3.145 ± 0.045 


3.498 ±0.048 


8.0 


0.703 ±0.035 


1.202 ± 0.040 


1.373 ± 0.046 


1.863 ±0.064 


2.099 ± 0.099 


3.047 ± 0.072 


3.375 ±0.073 



V 
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- no offsets 
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Figure 6. The median colours of ETGs, computed at R = R e , 
are plotted in the form of (g — X) — (g — X) (with X = rizY H K) 
as a function of the typical wavelength of filter X, where (g — X) 
is the median of all g — X values at different galaxy radii and 
masses (Tables [T] [2] and [3} . Different symbols refer to galaxies 
with different mass, as indicated in the lower-right of the plot. 
Error bars are la uncertainties on median colours. Dotted curves 
are the best-fitting SSP models from CB*MC (see text). Notice 
the discrepancy between models and data. The dashed curves are 
best-fitting SSP models after suitable correction offsets (OF Fx) 
were applied, improving the matching between models and data. 
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and models in preparation by ICharlot fc Bruzuall l|2013l . 
hereafter CB*). The CB07 models use the same sets of 
stellar tracks and spectral libraries as the BC03 models 
(see BC03 for details), except for the thermally pulsing 
asymptotic giant branch (TP-AGB) stars. For these stars, 
C B07 follow the sem i -empi rical evolutionary prescri ptions 
bv iMarieo fc Girardil l|2007n and lMarigo et al.l (|2008l ). The 
CB* population synthesis models used in this paper are 
based on the stellar evolution models comp uted with up- 
dated input physics for stars up to 15M© bv lBertelli et alJ 
(2008). These tracks are available for metallicities Z = 
{0.0001, 0.0004, 0.001, 0.002, 0.004, 0.008, 0.017, 0.04, 0.07}, 
with the solar metallicity being Zq = 0.017. In CB* the 
evolution of TP-AGB stars follows a recent prescription 
by Marigo & Girardi (private communication), which has 
been calibrated using observations of AGB stars in the 
Magellanic Clouds and nearby galaxies (|Girardi et al.ll2010l ; 
Mel bourne et al.ll2012h . In the optical range, the CB* mod- 



Figure 7. The colours used in this paper (g — X) are shown as 
a function of the effective wavelength of filter X. Top, middle, 
and bottom panels correspond to low-, intermediate-, and high- 
mass galaxies, respectively. For each mass bin, different galac- 
tocentric distances, R/R e = {0.1,0.5,1,2,4,8}, are plotted with 
different colours and symbols (see lower- left of top panel). The 
curves are best-fitting SSP colours for CB*MC models. Notice 
the good matching of observed and model colours. 



els are available for the I ndoUS llValdes et alJ 120041). Miles 
JSanchez-Blazquez et all I2006J), Stelib |Le Borgne et al.l 
l2003h . and BaSeL 3.1(| Westera et al.ll2002l ) spectral libraries. 
The NIR spectra of TP-AGB stars in CB * are s elected from 
the co mpilation bv lLancon fc Mouhcind (|2002l ), the IRTF 
library (|Ravner et all 20091 ). and the C-star model atlas by 
lAringer et all (|2009l ). For stars of metallicity near Zq , the 
number of available stellar spectra in the IndoUS and Miles 
libraries increases by a factor of roughly 20 with respect to 
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the Stelib library used in BC03. Thus, both in metallicity 
and spectral coverage, the CB* models represent a major 
improvement over the BC03 models, which translate into a 
better modelling of the NIR and optical galaxy colours. 

Recently, it has become clear that the treatment of the 
TP-AGB stars in the CB07 models tend t o overestimate the 
contribution by these stars in the NIR (|Kriek et al.l l20ld : 
iMelbourne et al.ll2012l ; IZibetti et ai1l2012l ). We use the CB07 
models in this paper for the sak e of completeness and com- 
paris on with previous work (e.g. lLa Barbera fc de Carvalhq 
12009] , and Paper IV). 

The simple stellar population (SSP) models adopted 
in the present study, and the labels we use to refer to 
them, are summarized in Tab. [4] The IMFs used are Scalo 
(BC03), Chabrier (CB07), while we consider both Chabrier 
and SalpeterQ (1955) IMF, for CB*. Models in the metal- 
licity range between 1/50 and 2.5Z© are considered, in the 
age range from 1 to 18 Gyr. The lower age limit of 1 Gyr 
is because we study ETGs, dominated by old stellar popu- 
lation, while the upper limit of 18 Gyr is the maximum age 
value available in the models 0. 

For a given model, synthetic colours are computed for 
three different star formation histories: 

- Simple Stellar Population. Equivalent to an instanta- 
neous burst of star formation. 

- t models. We assume an exponentially declining star for- 
mation rate with e-folding time r. The age of the model (i.e. 
the first epoch of star formation) and its metallicity are in 
the same ranges as for the SSPs. The r is varied from zero 
to the age of the model. 

- Burst model. The star formation rate is constant, with 
duration At. The age and metallicity of the model are varied 
in the same ranges as for r models. At is varied from zero 
to the age of the model. 

To compute synthetic colours, (g — X)mod {X — rizYHK, 
see below), each SP model is converted to the median red- 
shift of ETGs in our sample (z = 0.0614), and integrated 
with the grizYHK throughput curves (see Paper I for de- 
tails) . 

5.2 Colour fitting 

For a given galaxy mass bin, and given galactocentric dis- 
tance, we fit the corresponding median colours of ETGs by 
minimizing the following expression: 

2 _ Kg - X )o b s ~ (9 - X)mod + OFFxf 

x - 2^ a 2 w 

x 3-x 

where the summation extends over the available filters 
(X = rizYHK) , and the subscripts OBS and MOD re- 
fer to observed (Tables [T] [2] and [3]) and synthetic (model) 



i.e. a power-law IMF with slope 1.35, and lower and upper mass 
cutoffs of 0.1 and 100 Mq, respectively. 

2 We do not set the age of the Universe as an upper limit for the 
present analysis, as the absolute matching of model and observa- 
tions is still very uncertain. In fact, large differences exist between 
best-fitting ages from different models (see Sec. 16. ill , while rela- 
tive differences in age and metallicity (i.e. the radial gradients) 
are much better constrained. 



colours (see Sec. 15. ip . respectively. The terms OF Fx (X = 
rizYHK) are photometric zero-point offsets, adopted to im- 
prove the data and model matching (see below). In practice, 
the minimization consists of searching for the set of stellar 
population parameters (e.g. the age and metallicity for SSP 
models) that minimize Eq. [3] We remind the reader that 
g — J colours are not used in the minimization procedure, 
for the reasons explained in Sec. [4] Because some model 
parameters are intrinsically degenerate when fitting broad- 
band colours (e.g. the formation epoch and r, or At), for 
each model we give only the mass-weighted age parameter 
(Age, hereafter simply referred to as the age of the model), 
and metallicity (Z) and analyze how these quantities change 
with galaxy mass and galactocentric distance. This allows us 
to verify how the assumption of a given SFR affects the es- 
timate of age and metallicity profiles in ETGs. 

The terms OF Fx in Eq. [3] depend only on the adopted 
SP synthesis code (although the dependence is small, see 
below), and account for the fact that, in general, model 
colours do not match exactly the observed colours, because 
of, e.g., uncertainties in the photometric calibration and/or 
SP synthesis ingredients. To illustrate this point, we consider 
CB*MC models and plot in Fig. [6] the R = 1 R e colours of 
low-, intermediate-, and high-mass ETGs. For a clearer fig- 
ure including all colours in the same plot, we have subtracted 
the median colours, (g — X), the medians being computed 
among all mass bins and galactocentric distances. The val- 
ues of both (g — X) and OF Fx, for SSP models from differ- 
ent synthesis codes, are reported in Tab. [6] Dotted curves 
in Fig. [S] are obtained for CB07SC SSP models, by mini- 
mizing Eq. [3] with the values of OFFx set to zero. Notice 
the significant difference between observed and (best-fitting) 
model colours, ranging from a few hundredths of magnitude 
(g — r and g — K) up to ~ 0.1 mag in g — Y . To derive 
the offsets, OFF x , we set OFF x = in Eq. [3] and fit 
the observed colours with a given set of models (e.g. SSPs), 
for each galaxy mass bin and galactocentric distance. The 
values of OFFx are defined as the median differences of all 
observed and best-fitting colours. After including the OFFx 
terms in Eq. |3j and repeating the minimization procedure, 
we find that the median differences of observed and best- 
fitting colours become smaller than ~ 0.007 mag, proving 
that the entire procedure is self-consistent. We note that by 
using r and burst (rather than simple) SP models, the off- 
sets change less than ~ 0.01 mag, while the offsets differ by 
;$ 0.03 mag for SSPs among different SP codes (see Tab. [6]). 
The dashed curves in Fig. [6] show the offset-corrected model 
colours, i.e. (g — X)mod — OFFx, hereafter simply referred 
to as best-fitting model colours. Notice how the matching 
of observed and model colours improves significantly when 
applying the offsets. Fig. plots, as an example, the ob- 
served colours at all fiducial galactocentric distances (see 
Sec. [3}, for low-, intermediate-, and high- mass galaxies, and 
corresponding best-fitting colours for CB*MC SSPs. Over- 
all, we find that all the adopted models (i.e. population syn- 
thesis codes, and star formation histories) provide a fairly 
good match to the observations in all available passbands 
(but for J band, see above). We characterize the quality 
of colour fitting by computing the median rms value of dif- 
ferences between best-fit and observed colours, the median 
being computed among all fiducial radii. These rms values 
are summarized, for all colour fits performed in the present 
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Table 4. Stellar population synthesis models adopted in the present study. 



model stellar library IMF 



BC03 


STELIB 


Scalo 


CB07SC 


STELIB 


Chabrier 


CB*IC(CB*IS) 


INDOUX 


Chabrier (Salpeter) 


CB*MC(CB*MS) 


MILES 


Chabrier (Salpeter) 


CB*SC(CB*SS) 


STELIB 


Chabrier (Salpeter) 


CB*BC(CB*BS) 


BaSel 3.1 


Chabrier (Salpeter) 



work, in Tab. [5] Notice that the typical rms amounts to 
~ 0.02 mag. For what concerns best-fitting age and metal- 
licity, averaging over all radial bins and all models, the (ab- 
solute) values of age (metallicity) differ by ~ 18% (14%) 
when (i) one applies the OFFx's and (ii) the OFFx's are 
set to zero. On the other hand, relative values of age and 
metallicity (i.e. the gradients, see Sec. [6]) are virtually unaf- 
fected by the OFFx's. In order to estimate uncertainties on 
best-fitting age and metallicity, we perform 200 bootstrap 
iterations, where, each time, observed colours are shifted 
according to their errors and the entire minimization proce- 
dure is repeated. This bootstrap procedure provides a joint 
probability distribution function (PDF) for age and metallic- 
ity. Confidence intervals on age and metallicity are obtained 
at the 1 a equivalent, by computing the 16-th and 84-th 
percentiles of the marginalized PDF. 



6 AGE AND METALLICITY GRADIENTS 
OUT TO 8R B 

6.1 SSP models 

Fig. [8] plots the best-fitting age and metallicity as a func- 
tion of galactocentric distance, for BC03, CB07, and CB* 
(Chabrier IMF) SSP models. In the radial range from 
0.1 to lR e , the Age parameter (i.e. the "SSP-equivalent" 
age) increases, while metallicity decreases, implying that 



gradients, consistent with previous finding 


s (e.g. McChire 


19691; iPeletier et all 


19901; iGoreas et all 


19901; Gonzalez 


19931; iDavies et all 


1993; iFerreras et al 


l2005al. 2009; 


Sanchez-Blazauez et al.l 2006; La Barbera & de Carvalho 


20091; IClemens et all 


2009; Suh et al.l 1201 


x, iTortora et al. 


2010). In the present 


study, we supersede previous works, 



by exploring, for the first time, the behaviour of age and 
metallicity at large galactocentric distances. Overall, we find 
that at R > R e , Age keeps increasing, while metallicity de- 
creases, out to ~ 8 R c . Notice that the error bars on Age are 
on average large (up to 40% for some models and mass bins) 
at R > R c , reflecting the uncertainties on galaxy colours, 
as well as the fact that "old" (i.e. ^ 4 Gyr) ages are in- 
trinsically more difficult to constrain. This is due to the 
fact that for an old stellar population, a variation in age 
does not change much the model colors. Hence, the large 
uncertainties on Age at R > R c do not affect significantly 
those in metallicity, which remain reasonably small (< 10- 
20%) even at large galactocentric distances. Furthermore, 
the radial trend with metallicity appears quite robust across 
the different population synthesis models considered. Fig. [9] 
shows the age and metallicity profiles for CB* models when 



using a Salpeter IMF. Overall, the profiles are very similar 
to those presented in Fig. [8] implying that the choice of IMF 
in the stellar population models does not change at all our 
conclusions. 

In order to quantify the change in the slopes of the age 
and metallicity profiles, i.e. the age and metallicity gradients 
(hereafter Vt and V z), as a function of radius, we perform 
an orthogonal linear regression of the profiles in the inner 
(log(R/R c ) ^ 0) and outer (log(R/R c ) > 0) regions sepa- 
rately. The slopes give the inner and outer metallicity gra- 
dients (Vz,i and Vz,o), and the inner age gradient (Vt,i). 
Since age profiles become noisier at R > R e , we do not per- 
form any fit of the Age profile in the outskirts. Instead, we 
compute the average, logarithmic difference in Age, At, be- 
tween the outer region □ (R > 1R C ) and the galaxy centre 
(R = 0.1 Rc). On the contrary, errors on metallicity remain 
small even at large galactocentric distances (see above), al- 
lowing us to estimate meaningful outer metallicity gradi- 
ents. Errors on Vt, Vz, and At are estimated as the errors 
on age and metallicity (see end of Sec. 15. 2[) , performing 200 
bootstrap iterations, where colours are shifted according to 
their uncertainties, and age and metallicity gradients are re- 
estimated, by repeating the \ 2 minimization procedure at 
each fiducal radius. For all different SSP models, the val- 
ues of Vz,i and Vz,o (Vt,i and At) are reported in Tab. [7] 
(Tab. HO). 

Fig. [8] unsurprisingly shows that absolute estimates of 
Age and metallicity differ significantly among different mod- 
els, while relative trends, i.e. with respect to mass and galac- 
tocentric distance, do not (but see below). In the galaxy 
centre (R ~ 0.1 R e ), BC03 models predict a super-solar 
metallicity (~ IAZq) and an Age of ~ 5 Gyr, while for 
CB07, the central metallicity is sub-solar (~ O.65Z0), and 
Age older (with respect to BC03) by ~ 4 Gyr. The new 
CB* models also exhibit significant reciprocal differences in 
Age (up to a few Gyrs), and relatively small differences in Z 
(~ 0.1 dex), when comparing results based on different stel- 
lar libraries (CB*IC; CB*MC; CB*SC; CB*BC). Notice that, 
for high-mass ETGs, all CB* models predict either solar 
(CB*MC; CB*BC) or super-solar (CB*IC; CB*SC) metallici- 
ties. Hence, in contrast to the CB07 models, the photometric 
constraints of the CB* models are more consistent with the 
results fro m spectroscopic studies of the central reg ions of 
ETGs (e.g. iGallazzi et al.ll2006l ; iPasauali et al.ll2010f i. Since 
the absolute Age is strongly model dependent (even when 
changing only the adopted stellar library in CB*), we have 



3 In practice, we compute the mean value of Age among the 
three fiducial radii of 2, 4, and 8 R c , and then take the logarithmic 
difference between this mean value and the Age at R = 0.1 R e . 
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Table 5. Median rms of differences between best-fit and observed colours, for all models and samples analyzed in this work. The median 
values are computed among all fiducial radii, and inform on the quality of the colour fitting. For CB* models, the offsets corresponding to 
a Salpeter (S), rather than Chabrier (C), IMF are given in brackets. 



model 




median rms 






low 


intermediate- 


high.-m.ciss 


BC03 SSPs 


0.017 


0.015 


0.008 


CB07SC SSPs 


0.016 


0.020 


0.013 


CB*IC (CB*IS) SSPs 


0.022(0.020) 


0.017(0.012) 


0.019(0.019) 


CB*MC (CB*MS) SSPs 


0.032(0.016) 


0.010(0.010) 


0.015(0.015) 


CB*SC (CB*SS) SSPs 


0.015(0.015) 


0.014(0.015) 


0.014(0.012) 


CB*BC (CB*BS) SSPs 


0.019(0.016) 


0.013(0.011) 


0.024(0.022) 


CB*SC r models 


0.016 


0.015 


0.014 


CB*SC burst models 


0.016 


0.015 


0.013 


CB07SC (E(B -V) = 0.04) 


0.026 


0.024 


0.021 


CB07SC {E(B -V) = 0.08) 


0.023 


0.019 


0.020 


CB07SC (for field ETGs) 


0.020 


0.019 


0.026 


CB*MC (for field ETGs) 


0.015 


0.017 


0.024 


CB07SC (for group ETGs) 


0.019 


0.021 


0.020 


CB*MC (for group ETGs) 


0.012 


0.020 


0.017 



Table 6. Median colours of ETGs and photometric offsets, OF Fx, applied to different SSP models. For CB* models, the offsets corre- 
sponding to a Salpeter (S), rather than Chabrier (C), IMF are given in brackets. 



colour 


median 








OFF x 












BC03 


CB07SC 


CB*IC(S) 


CB*MC(S) 


CB*SC(S) 


CB*BC(S) 


9- 


- r 


0.803 


-0.033 


-0.019 


-0.028(-0.027) 


-0.022(-0.022) 


-0.025(-0.022) 


-0.016(-0.018) 


9~ 


- i 


1.197 


0.032 


0.031 


0.028(0.029) 


0.035(0.034) 


0.024(0.023) 


0.039(0.039) 


9~ 


- z 


1.458 


0.044 


0.016 


0.054(0.052) 


0.042(0.034) 


0.049(0.044) 


0.023(0.018) 


9- 


Y 


2.152 


-0.119 


-0.131 


-0.109(-0.114) 


-0.129(-0.132) 


-0.117(-0.121) 


-0.131(-0.131) 


9 - 


H 


3.286 


0.035 


0.060 


0.016(0.018) 


0.033(0.040) 


0.028(0.032) 


0.044(0.045) 


9 - 


K 


3.627 


0.037 


0.041 


0.014(-0.011) 


0.017(0.027) 


0.017(0.024) 


0.025(0.034) 



not set the age of the Universe as an upper limit for the mod- 
els (see Sec. l5.1[l . and several profiles in Fig.[8]keep increasing 
outwards up to the maximum allowed value of 18 Gyr. Im- 
posing an Age upper limit of 13.75 Gy r (i.e. the current best 
estimate for the age of the Universe, iKomatsu et all 1201 lh 
would produce an upper cutoff in the Age profiles, while not 
affecting significantly the metallicity profiles. Considering 
the large Age differences among different models, enforcing 
an upper limit on Age is not meaningful for the models, so 
we decided not to apply this constraint. 

Overall, the key features of the age and metallicity 
trends are summarised as follows: 

- Metallicity is a monotonic, decreasing function of galac- 
tocentric distance in ETGs, out to (at least) 8R e . This is 
the most robust result, holding true for all mass bins, and 
for all different models probed here. 

- The values of At in Tab. [8] are always positive. In most 
cases (and in particular for all CB* models with Chabrier 
IMF), this result is found above the 2<r level. For high- mass 
galaxies, the At is positive at more than the 3cr level for 
all models. We conclude that ETGs, in particular the most 
massive ones, host significantly older stellar populations in 
their halos (R ^, few x R e ) than in their cores. 

- In the inner region, the metallicity gradient, Vz,i, steep- 
ens at high mass, while Vt,i increases with M*. These trends 
hold for all models, except for CB*/IndoUS, where no mass 



dependence is detected. A steepening of Vz,i and Vt,i with 
galaxy mass has also been found in previous papers based o n 
the SPIDER survey (e.g. Paper IV: lLa Barbera et al-lEm! ) . 
where colour gradients were fitted (rather than colours) with 
SP models. Notice that at high mass, the Vt,; is always pos- 
itive, regardless of the SP model, consistent with our find- 
ings for At in the external regions. At intermediate- and 
low-mass, the Vt is either consistent with zero, within the 
errors, or significantly positive. In particular, at low mass, 
the age trend is flat for most models (i.e. BC03, CB07SC, 
CB*SC, and CB*BC). 

- For low-mass ETGs, the metallicity gradient is steeper 
at large radii than in the centre, i.e. Vz,o<Vz,i. While 
this result holds for all models, it is statistically signifi- 
cant (at more than a 2-3cr level) only for some models 
(e.g.CB*SC). At high-mass, the same steepening is observed 
only for some models when using a Chabrier IMF (e.g. 
CB*IC), and for all models with a Salpeter IMF (in which 
case, however, the steepening is statistically significant only 
for CB*IS models). If a bottom-heavy, Salpeter-like, IMF 
is m ore appropriate for high-, rel a tive to low-mass ETG s 
(e.g. Ivan Dokkum fc Conrov|[20Tl] : ICappellari et al ] |2012T I. 
we might be tempted to conclude that the metallicity gra- 
dient of ETGs steepens at large galacto-centric distances, 
independent of galaxy stellar mass. 
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Figure 8. Age and metallicity profiles obtained by fitting galaxy colours with SSP models, as a function of galactocentric distance, for 
low-(blue), intermediate- (green), and high-(red) mass ETGs. Different symbols are used for different galaxy mass bins, as illustrated 
in the top-left panel. From left to right, top and bottom panels refer to different models (BC03, CB07SC, CB*IC, CB*MC, CB*SC, 
CB*BC), as illustrated by the black labels in the top panels. All models correspond to a Chabrier IMF. Error bars denote the 16th and 
84th percentile uncertainties (corresponding to la errors for a normal deviate) on best-fitting age and metallicity. 




Figure 9. Same as Fig.[8]but plotting only the results for the CB* models using a Salpeter IMF. 



6.2 Exponential and burst models 

Fig. [10] shows radial trends of (mass-weighted) age and 
metallicity as obtained by fitting galaxy colours with com- 
posite stellar populations (r and burst models). For brevity 
reasons, we only show one illustrative case, corresponding to 
CB*SC models (Sec. 15.21) . Both Age and metallicity exhibit 
similar trends to those obtained with SSPs (Fig. [8}. In par- 
ticular, we see that (i) Age increases outwards, being oldest 
in the galaxy halo region; (ii) the age gradient is stronger 
(more positive) at higher galaxy mass , with low-mass ETGs 
having a flat age trend out to ~2-4R c ; and (iii) metallicity 
decreases outwards, with a steeper (more negative) gradient 
(especially for low- and intermediate-mass galaxies) at large 
galactocentric distances. One may also notice that individ- 
ual values of Age are somewhat more sensitive (in contrast 



to metallicity) to the adopted star formation history (i.e. 
SSP, t, and burst models). For instance, at R > 3R C , and 
for low- mass ETGs, both burst and SSP models (see Fig. [8| 
give the oldest available Age (18 Gyr) in the input grid 
of values, while for r models the Age is ~ 13 Gyr (see blue 
curves in the upper panel of the Figure). On the other hand, 
all models agree very well (within error bars) on the inferred 
value of metallicity at each radial position. 



6.3 Internal reddening 

So far, we have neglected the role of internal reddening, 
which mi ght mimic the presence of metallicity gradients 
in ETGs iGoudfrooii fc de Jond 1 19951 ; ISilva fc Wise] Il996l ; 
I Wise fc Silvalll996l ). Indeed, if dust gradients are present, 
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Figure 10. The best-fitting age (top) and metallicity (bottom) 
are shown as a function of galactocentric distance, as in Fig. \E\ 
for low- (blue), intermediate- (green) , and high- (red) mass ETGs. 
Smaller (larger) symbols and solid (dashed) curves refer to CB*SC 
exponentially declining (burst) models. Notice that the trends are 
very similar to those obtained with SSP models (Fig. [HJ. 



Figure 11. The same as right panel of Fig . [8] but comparing 
results of CB07SC with ICardelli et aLl j 19831 extinction law ap- 
plied, with a colour excess of E(B — V) = 0.04 (solid curves) and 
E(B — V) = 0.08 (dashed curves). Notice that internal reddening 
does not affect significantly the trends of age and metallicity with 
galactocentric distance. 



a correlation of colour gradients and the amount of inter- 
nal extinction should be detected. Such correlation does not 
exist but for a small fraction (< 10%) of ETGs in our sam- 
ple (see Paper IV). Also, previous studies have found nei- 
ther a significant correlation be tween colour gradients and 
IRAS 100 n flux (|Michardl20"05h . n or a contribution of dif- 
fuse dust to the colour gradients l|Savov et all [2009). Al- 
though the evidence supports our assumption of neglecting 
extinction gradients in ETGs, the presence of uniform in- 
ternal reddening inside the galaxies may still affect our con- 
clusions, as it could bias the best-fitting values of age and 
metallicity at each radial position. To address this issue, we 
have applied dust extinction corrections to o ne of the avail- 
able s ets of CB* SSPs (CB*MC), adopting the lCardelli et all 
( 1989) extinction law, and repeated the SSP fitting proce- 
dure for two cases, i.e. a colour excess of E(B — V) = 0.04 
and E(B — V) = 0.08. These values correspond to the me- 
dian and the 75% level of the probability distribution of the 
colour excess parameter, i nferred by running the sp ectral 
fitting code STARLIGHT (jCid Fernandes et alJliooa ). with 
the CB07SC models, on the SDSS spectra of our sample of 
ETGs (see Paper V for details). Fig. 1111 shows the age and 
metallicity profiles obtained when adopting the above red- 
dening values. The trends of age and metallicity are fully 
consistent with those obtained for the dustless models (see 
CBTVIC profiles in Fig. [8]). 



6.4 Comparison with recent findings 

O ur findings (e.g. Fig. [St can be com pared with those 
of ICoccato. Gerhard fc Amaboldil (|2010l . hereafter CGA10), 
who derived radial profiles of stellar population properties 
(i.e. SSP-equivalent age, [Z/H], and [a/Fe]), based on spec- 
troscopic data, for one of the two brightest galaxies in the 
Coma cluster (NGC 4889), out to a galactocentric distance 
of ~ 4R C . CGA10 found that the stellar populations in the 
outer halo of NGC 4889 tend to be older than those in the 
inner regions, in agreement with our results. At R ^ 1.2 R c , 
the metallicity profile of NGC 4889 flattens, while [a/Fe] 
shows a steep negative gradient. Although the optical+NIR 
photometry allows us to reasonably disentangle the average 
contribution of age and metallicity to galaxy colours (see 
Paper IV), we are not able to constrain both [Z/H] and 
[a/Fe] separately, and our metallicity profiles can indeed be 
affected by radial changes of [a/Fe]. 

For instance, the steepening of the metallicity profile 
at R > R c in Figs. [5] and [TU] may be due to a steepening 
of the [a/Fe] profile at large galactocentric distances, as ob- 
served by CGA10, or, alternatively, it can be a true change in 
[Z/H], implying that the result of CGA10 does not apply to 
the entire population of ETGs in our sample. The fact that 
we obtain similar trends of Age and [Z/H] for both r and 
SSP models favours a true variation of the metallicity profile 
in the galaxy outskirts rather than a variation of the star 
formation time-scale (i.e. abundance ratios) at large radii. 
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Table 7. Mctallicity gradients (V z) °f ETGs estimated using different SSP models. Subscript i and o refer to the inner (R R e ) and 
outer (R e < R ^ 8R C ) radial ranges where the gradients are estimated. 



model METALLICITY GRADIENTS 













V z , 






low- 


intermediate- 


high-mass 


low- 


intermediate- 


high-mass 


BC03 


-0.31 ±0.06 


-0.52 ±0.05 


-0.55 ± 0.04 


-0.38 ±0.17 


-0.62 ± 0.12 


-0.29 ±0.12 


CB07SC 


-0.32 ±0.04 


-0.38 ±0.03 


-0.41 ± 0.03 


-0.57 ±0.26 


-0.80 ±0.13 


-0.74 ±0.20 


CB*IC 


-0.40 ± 0.04 


-0.44 ±0.05 


-0.39 ± 0.05 


-0.46 ±0.06 


-0.54 ±0.07 


-0.55 ±0.12 


CB*MC 


-0.30 ±0.04 


-0.35 ±0.05 


-0.49 ± 0.07 


-0.58 ±0.14 


-0.43 ± 0.10 


-0.32 ±0.13 


CB*SC 


-0.29 ±0.05 


-0.39 ±0.07 


-0.46 ± 0.05 


-0.57 ±0.09 


-0.50 ±0.12 


-0.54 ±0.15 


CB*BC 


-0.27 ±0.03 


-0.34 ±0.06 


-0.41 ± 0.07 


-0.61 ±0.12 


-0.54 ±0.11 


-0.51 ±0.13 


CB*IS 


-0.39 ±0.03 


-0.43 ±0.04 


-0.38 ± 0.06 


-0.49 ±0.07 


-0.55 ± 0.13 


-0.62 ±0.08 


CB*MS 


-0.30 ±0.03 


-0.33 ±0.03 


-0.41 ± 0.08 


-0.64 ±0.13 


-0.58 ± 0.13 


-0.51 ±0.13 


CB*SS 


-0.29 ±0.05 


-0.37 ±0.05 


-0.43 ± 0.06 


-0.65 ±0.10 


-0.54 ±0.13 


-0.58 ±0.16 


CB*BS 


-0.29 ±0.03 


-0.33 ±0.03 


-0.39 ± 0.06 


-0.59 ±0.11 


-0.63 ± 0.20 


-0.59 ±0.16 



Table 8. Same as Tab. [7] but for age, rather than metallicity, gradients (Vt). Notice that instead of computing age gradients in the 
galaxy outskirts, we measure the mean logarithmic age difference, At, between the outer radial range and the galaxy centre (see text for 
details). 



model AGE GRADIENTS 



Vt,» At 





low- 


intermediate- 


high-mass 


low- 


intcrmcdiatc- 


high-mass 


BC03 


0.01 ± 0.03 


0.19 ±0.03 


0.26 ±0.02 


0.13 ± 0.06 


0.31 ±0.05 


0.26 ± 0.04 


CB07SC 


0.03 ± 0.03 


0.09 ±0.02 


0.15 ±0.02 


0.03 ± 0.08 


0.25 ±0.05 


0.24 ±0.07 


CB*IC 


0.16 ± 0.03 


0.22 ±0.04 


0.18 ±0.06 


0.25 ± 0.08 


0.42 ± 0.08 


0.26 ± 0.10 


CB*MC 


0.04 ± 0.04 


0.10 ±0.05 


0.34 ±0.05 


0.28 ± 0.07 


0.19 ±0.09 


0.47 ±0.08 


CB*SC 


-0.00 ±0.05 


0.12 ±0.06 


0.26 ±0.05 


0.27 ±0.09 


0.20 ±0.09 


0.32 ± 0.09 


CB*BC 


-0.04 ±0.05 


0.09 ±0.05 


0.19 ±0.06 


0.18 ± 0.08 


0.33 ±0.08 


0.36 ± 0.10 


CB1S 


0.15 ± 0.05 


0.21 ±0.05 


0.17 ±0.09 


0.25 ± 0.05 


0.43 ± 0.07 


0.37 ±0.10 


CB*MS 


0.04 ±0.06 


0.08 ±0.06 


0.20 ±0.09 


0.27 ±0.06 


0.32 ±0.09 


0.38 ± 0.10 


CB*SS 


-0.00 ±0.06 


0.09 ±0.08 


0.21 ±0.07 


0.28 ± 0.06 


0.20 ±0.08 


0.33 ± 0.09 


CB*BS 


-0.01 ± 0.06 


0.03 ±0.05 


0.16 ±0.05 


0.11 ±0.09 


0.26 ±0.07 


0.34 ± 0.07 



In this regard, our re sults are mor e cons istent with those 
recently presented by I Greene et al. (2012), who found, for 
a sample of eight massive ETGs, that metallicity gradients 
seen within the effective radius continue to fall to outer radii 
of ~ 2.5 R , while in constrast, [a/Fe] does not drop substan- 
tially at large radii in these objects. Notice that so far, no 
significant a bundance ratio gradien t s have been detecte d in 
ETGs (e.g., lKuntschner et alJkoid : ISpolaor et alj|20ich . 



7 THE EFFECT OF GALAXY ENVIRONMENT 

Fig-[H]compares the radial profiles of age and metallicity for 
the subsamples of field and group ETGs (Sec. [2]). For each 
subsample, we derive the median (g — X) colour profiles 
(X — rizYHK) and fit them, at the fiducial galactocentric 
distances, as described in Sec. 15.21 We show here only the 
results of fitting C B07SC SSPs, for com parison with our 
previous work (e.g. lLa Barbera et al.ll201ll . hereafter LF11), 
and those obtained for one of the CB* models (CB*MC), 
as other models arrive at similar conclusions (see Sec. 16. ip . 
Solid and dashed curves in the Figure refer to field and group 
galaxies, respectively. Black symbols and error bars, on the 
left and right sides of each panel, represent the mean values 
of age and metallicity along with their error bars at the 
extreme fiducial points of 0.1 and 8R C , respectively. The 



mean values are obtained by averaging the values of age and 
metallicity for the three mass bins at the extreme fiducial 
points of 0.1 and 8 R e , respectively. For each subsample, we 
derive the inner and outer age (metallicity) gradients, Vt,j 
and At (Vz,i and Vz, ), as described in Sec. [6] Their values 
are reported in Tab. Hand Tab. [TO] for CB07SC and CB*MC 
models, respectively. 

Central ages. In the inner region (R 5j R e ), group ETGs 
have positive age gradients (age increasing outwards) for 
each mass bin (Vt,i> 0), with the strongest gradients at 
high mass. On the other hand, field ETGs exhibit flatter 
V t ,i's. For CB07SC SSPs, field galaxies vary from (slightly 
negative but consistent with) zero gradients at low-mass 
(Vt,i= —0.01 ± 0.06) to slightly positive at high mass 
(Vt,i~ 0.05), while for CBTVIC, the Vt.i is positive, but 
consistent with zero (at the ~ 1 a level, see Tab. I10[) for all 
mass bins. These findings are consistent with our previous 
results (LF11), based on a larger sample of SPIDER ETGs, 
where field ETGs (at low- and intermediate-mass) are found 
with negative - but sma ll - age gradients (consistent with, 
e.g., iFerreras et aT1l2009l ). 

Central metalhcities. No significant difference is detected, 
within the errors, between the Vz,i of field and group ETGs, 
although for CB*MC models, at high-mass, the metallicity 
gradient is significantly steeper for group (Vz,,~ —0.58) rel- 
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Table 9. Metallicity (Vz) and age (Vt) gradients of ETGs as a function of environment, for CB07SC SSP models. The subscripts i and 
o refer to the inner (R R e ) and outer (R c < R < 8R C ) radial ranges where gradients are estimated. 



mass range 


Vz,i Vz, 


V M 


A t 


3.00 < M*[10 10 M Q ] < 5.29 
5.29 < M*[10 10 A'/ Q ] < 9.42 
9.42 < Af*[lO lo M ] < 70.71 


FIELD ETGs 
-0.32 ±0.04 -0.87 ±0.25 
-0.35 ±0.04 -0.13 ±0.17 
-0.39 ±0.03 -0.15 ±0.22 


-0.01 ±0.06 
0.05 ± 0.05 
0.06 ±0.05 


-0.02 ±0.11 
0.01 ±0.11 
0.06 ±0.10 


3.00 < M*[10 10 M Q ] < 5.29 
5.29 < A/*[10 10 Af Q ] < 9.42 
9.42 < Af*[lQ 10 M Q ] < 70.71 


GROUP ETGs 
-0.36 ±0.04 -0.47 ±0.18 
-0.36 ±0.03 -0.92 ±0.12 
-0.41 ±0.03 -0.94 ±0.25 


0.08 ± 0.05 
0.10 ±0.04 
0.20 ±0.04 


0.10 ±0.08 
0.25 ±0.04 
0.30 ±0.10 



Table 10. The same as Tab.[9]but for CB*MC rather than CB07SC SSP models. 



mass range \7 z ,i Vz,o V tji At 

FIELD ETGs 

3.00 < M*[1O 1O M ] < 5.29 -0.37 ±0.06 -0.87 ±0.30 0.09 ± 0.08 0.38 ± 0.09 

5.29 < A/*[1O 1O M ] < 9.42 -0.34 ±0.08 -0.36 ±0.16 0.12 ±0.13 0.21 ±0.13 

9.42 < Af*[10 10 M Q ] < 70.71 -0.37 ±0.07 -0.15 ±0.16 0.11 ±0.08 0.16 ±0.06 

GROUP ETGs 

3.00 < A4*[1O 1O M ] < 5.29 -0.39 ± 0.06 -0.57 ±0.11 0.22 ±0.10 0.43 ±0.10 

5.29 < M*[1O 1O M ] < 9.42 -0.32 ± 0.03 -0.56 ±0.10 0.11 ±0.03 0.19 ±0.03 

9.42 < M*[iO lo M ] < 70.71 -0.58 ± 0.04 -0.52 ±0.13 0.52 ± 0.06 0.54 ± 0.06 



ative to field (Vz,,~ -0.37) ETGs. Note that LF11, based 
on larger samples of ETGs (i.e. smaller error bars on V z,i), 
found the inner metallicity gradient of group ETGs to be 
steeper (at all masses), by about —0.05, than that of field 
galaxies. On average, at 0.1 R e , group ETGs are more metal- 
rich than their field counterparts, with <5(log Z/Zq) — 0.05± 
0.02 (CB07SC) and S(\ogZ/Z Q ) = 0.08±0.03 (CB*MC), av- 
eraging over all three mass bins (see star symbols in the left 
side of each bottom panel in Fig. I12p . 

Outer ages. In the galaxy outskirts (R > R e ), for group 
ETGs, the Age keeps increasing outwards, with very old 
ages at the furthest radial distances probed. Averaging 
over all three mass bins, for group ETGs, we find (A t ) = 
0.25 ± 0.05 dex (CB07SC) and (A t ) = 0.43 ± 0.09 dex 
(CB*MC). On the contrary, for field ETGs, the Age trend 
is strongly dependent on the stellar population model one 
adopts to fit the color profiles. For CB07SC, Age does 
not increase significantly with galactocentric distance, with 
(At) = 0.01 ± 0.07 dex, while for CB*MC, Age increases 
outwards, with (At) = 0.21 ± 0.08 dex. Hence, field ETGs 
have either positive or null age gradients in the outskirts 
(depending on the model), while group ETGs have positive 
age gradients in the outskirts. 

Outer metallicities. For group ETGs, the metallicity gra- 
dient is found to steepen in the galaxy outskirts when us- 
ing CB07SC models ((Vz, D ) = -0.8 ± 0.1, compared to 
(Vzj = —0.38 ± 0.02, averaging all three mass bins). 
This result is somewhat model dependent, being confirmed 
- for CB* models - only at low- and intermediate-mass 
«Vz, ) = -0.57±0.07 compared to {V z,i) = -0.36±0.03). 
For field ETGs, no steepening is observed at intermediate- 
and high-mass, for both CB07SC and CB*MC models. At 
low-mass, galaxies in the field also exhibit a steeper gra- 
dient in the outskirts, for all models (Vz, ~ —0.87 with 



respect to Vz,;~ -0.32, for CB07SC, and V z , 4 ~ -0.37, for 
CB*MC). Hence, the steepening of the metallicity gradients 
in low-mass (M* ~ 4x 1O 1O M ) ETGs is detected in all envi- 
ronments, independent of the SP model one adopts to fit the 
color profiles (consistent with what found in Sec. l6.1|) . Notice 
that at 8 R c , group ETGs are on average less metal-rich than 
their field counterparts Q, in contrast to the observations in 
the inner region (see above). 



8 SUMMARY AND CONCLUSIONS 

We have used optical (SDSS-DR6) and near-infrared 
(UKIDSS-DR4) data for 674 massive ETGs to measure, for 
the first time, colour gradients out to ~ 8R e , and study 
how ETGs assemble their mass. The sample comprises sys- 
tems found in galaxy groups and in the field, with stellar 
masses ranging from ~3x 1O 1O M to ~ 7 x 1O 11 M . The 
critical systematic effects that may affect our conclusions 
have been carefully taken into account, and we show that 
none of them (e.g. sky subtraction, colour dependence of the 
PSF outer wings, and the approach adopted to derive the 
colour profiles) introduce any systematic biases. The stacked 
colour profiles show remarkable linearity out to R ~ 8R C , 
getting bluer in the outer regions, irrespective of the mass 
bin, with the exception of g — r where a curvature is ap- 
parent. Nevertheless, this effect is small and does not affect 



4 This is truly the case only for high- and intermediate-mass 
ETGs. In fact, as one can see in the bottom panels of Fig. 1121 the 
green and red dashed lines falls systematically below the green 
and red solid lines at R/R e > 2. At low-mass, group ETGs are 
more metal-rich than their field counterparts, at all masses. 
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Figure 12. Comparison of age (upper panels) and metallicity (lower panels) profiles of field (solid curves) and group (dashed curves) 
ETGs. Left and right panels refer to CB07SC and CB*MC models, respectively (see labels in the lower-right of top panels). For R = 0.1 R c 
(R = 8R C ), the mean values of age and metallicity are plotted with black star symbols on the left (right) side of each panel, with small 
and big stars corresponding to field and group galaxies, respectively. Black error bars are the mean errors on age and metallicity at 
R = 0.1 R c and R = 8R C . 



the stellar population analysis, as the deviations are com- 
parable to the random errors at the radii where the curva- 
ture is detected. For the most massive galaxies, the stacked 
g — r colour profile is c o mpar ed with that recently obtained 
bv lTal fc van Dokkuml l|201ll) . finding good agreement. 

Colour profiles are translated into SSP-equivalent Age 
and metallicity gradients using a variety of stellar population 
models, assuming either a Chabrier or Salpeter IMF, as de- 
tailed in Section 5.2. In particular, we exploit the new, state- 
of-the-art, stellar population models of Chariot & Bruzual 
(2013, in preparation), based on different stellar libraries 
(IndoUS, Miles, STELIB, BaSeL). Errors on both age and 
metallicity are carefully estimated by accounting for uncer- 
tainties on observed colors. We find very unambiguously that 
metallicity decreases from 0.1 out to 8R e (negative gradi- 
ent), regardless of the mass bin considered (see Tab. 6). 
This result is independent of the adopted stellar population 
model and IMF. In addition, we find that the metallicity 
gradient tends to steepen at large galacto-centric radii. The 
significance of this result changes with the adopted stellar 
population models and galaxy mass bin (being more sig- 
nificant at low mass). As for the Age parameter, the situ- 
ation is more complex. From 0.1 to lR e , we find that for 



intermediate- and high-mass ETGs, Age increases with ra- 
dius, while at low mass, no clear trend with galacto-centric 
dist ance is detected, consist ent with our previous results (see 
e.g.|l][ Barbera et al.ll2~01ll and references therein) . The new 
result of the present work is that in the outer regions, from 1 
to 8 R e , the stellar populations of ETGs are even older than 
in their centres. This is more significant for galaxies in the 
highest stellar mass bin (i.e. M* Si lO n M0). We emphasize 
that using models with extended star formation histories, 
like exponentially declining (r) and finite burst models, in- 
stead of SSPs, yields similar results within the quoted error 
bars. We also confirm that the presence of uniformly dis- 
tributed internal reddening (dust) does not affect the trends 
of age and metallicity with galactocentric distance. 

Finally, we analyze how the age and metallicity profiles 
depend on the environment where galaxies reside. Group 
ETGs have positive age gradients out to ~ 8R C , with a 
very old stellar population in the outskirts. Their metallic- 
ity gradients steepen at R > 1 — 2R e . In constrast, field 
ETGs do not have significant age gradients in the inner re- 
gion, while at R £ 1 R c , the Age shows either a flat trend or 
increases outwards (depending on the model). The metal- 
licity gradient of field ETGs is found to steepen signifi- 
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cantly in the outskirts only for low-mass systems. Notice 
that the results for group galaxies are essentially the same 
as those for the whole sample of ETGs (see above), consis- 
tent with the fact that most of our sample (~ 80%) com- 
prises ETGs either residing in (or close to) galaxy groups 
(see Sec. [2| . Our findings point to the importance of aper- 
ture effects when comparing the stellar population content 
of ETGs in different environments, as in the galaxy out- 
skirts, group galaxies are older (or coeval, depending on the 
model) and less metal-rich than field ETGs, while the op- 
posite holds in the central regions, where group ETGs are 
found to be more metal-rich than their field counterparts. 
This might explain, at least in part, why different results 
have been reported in the literature about the environmen- 
tal dependence of the metal content of ETGs (see Paper 
III), with di fferent studies r eport in g field ETGs to b e more 
metal-rich dThomas et al l 120051; Ide La Rosa et ail l2007t 



Kuntschner et all 12002]; IClemens et all 120091: IZhu et al 



201C ), as metal-rich as ~(|Bernardi et al.ll2006l " lAnnibali et al 



20071 ), or even more metal -poor than their cluster counter- 
parts (|Gallazzi et al.ll2006l ). 



We have shown here that the stellar populations in the 
halos of ETGs are consistently older and more metal poor 
than in their cores. This trend depends on mass and, to a 
lower degree, environment, but overall, the age and metal- 
licity radial gradients are similar for all ETGs in the sample 
(i.e. more massive than 3 x 10 10 Mq). Within the core of 
the galaxies, the observed metallicity gradients - between 
Vz,j= —0.3 and —0.4 - are consistent with the simula- 
tions of iKobavashil (|2004h for non-major merger systems. 
For galaxies where major-merging is significant, her simu- 
lations give shallower metallicity gradients, around —0.2. It 
is worth noting that our metallicity gradients get slightly 
steeper with increasing mass, supporting the scenario of 
fast and early formation for the most massive ETGs. Ex- 
tending the analysis radially out to 8R e allows us to ex- 
plore the different channels of star formation and assem- 
bly inFJTGs. Our r esults are consistent with the scenario 
of lOser et ail (|2010l ) , where two phases of formation op- 
erate: an early-phase of in-situ star formation - building 
the core of the galaxy - followed by the accretion o f smal l 
satellites. The recent simulations of lLackner et all (120121 ') 
quantify in more detail the buildup of the stellar compo- 
nent, concluding that major mergers do not dominate the 
accretion history of massive galaxies, with smaller systems 
consisting of older, and more metal-poor stellar popula- 
tions building up the outer halo, consistent with our anal- 
ysis. Their simulations suggest age (metallicity) differences 
of At ~ 2.5 Gyr (AlogZ ~ — 0.15dex) between accreted 
stars and those formed in-situ. This process is fully con- 
sistent with the observed evolution of massive early-type 
galaxies on the ma s s-size plane (see, e.g. iDaddi et al.ll2005l ; 



iTruiillo et all 120061 . I2OIII ; Ivan Dokkum et alj|2008h . where 
the core is quickly fo rmed during an early phase (see, e.g. 
iFerreras et alll2009bl) , fo llowed by the grow th of the halo via 
minor mergers (see, e.g. iNaab et al.l l2009). 



APPENDIX A: COMPARISON OF 
PARAMETRIC AND NON-PARAMETRIC 
LIGHT PROFILES 

Fig. IA1I compares parametric and non-parametric stacked 
light profiles for each of the three galaxy mass bins (from 
left to right), and each waveband, X — grizYHK (from 
top to bottom), used in the present study. We do not con- 
sider J band, as this is not used for the stellar popula- 
tion analysis (see Sec. 0] for details). For each galaxy im- 
age, and a given waveband X, we measure mean surface 
brightness values on a set of concentric circles, with radii 
equally spaced by 0.5 pixels, out to a maximum distance of 
~ 110 arcsec from the galaxy center. All objects around the 
given galaxy are masked out wi t h the software 2DPHOT, as 
described in lLa Barbera et ail (|2008l ) . The surface bright- 
ness profiles of all galaxies are scaled by the correspond- 
ing iie.x's, and normalized to have the same flux (within 
an aperture of 2R Bt x)- The so-normalized profiles are aver- 
aged within radial bins logarithmically spaced by 0.15 dex. 
This procedure provides non-parametric light profiles (not 
corrected for PSF), with grey regions in Fig. I All marking 
the 1 a mean confidence contours around them. The same 
procedure, applied to the PSF-convolved best-fitting two- 
dimensional Sersic models, gives the parametric light pro- 
files, plotted as black curves in Fig. I All Overall, we find a 
good agreement between observed and model profiles, out 
to the largest galactocentric distances considered in this 
work (~ 8 R e ). This supports previous claims that a sin- 
gle Sersic law is able to describe a ccurately the light dis- 
tribution of ETGs out t o ~ 8 R e l|Kormendv et al 1 120091 ; 
iTal fc van DokkumllioriT) . supporting the robustness of our 
parametric colour estimates over a wide radial range. 
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